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Abstract 

Six novel polymorphic microsatellite markers were isolated from enriched libraries in Myrica faya Ait., recently re- 
named Morella faya, (fayatree, firetree, or firebush) in order to examine the genetic diversity in natural populations. 
Also, test cross-specific amplification and genetic diversity in Myrica rivas-martinezii, which is endemic on the Ca- 
nary islands. Microsatellite loci were screened in 225 individuals of both species from different islands of the 
Canarian archipelago. All markers were successfully amplified from both Myrica species, with an average number of 
6.5 and 9.3 alleles per locus in M. rivas-martinezii and M. faya, respectively. There was no evidence for linkage dis- 
equilibrium between loci, and the probability of null alleles ranged from 0.01 to 0.17. 
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M. faya is native to the northern islands of Maca- 
ronesia, the Azores, Madeira, and the Canaries in the North 
Atlantic Ocean. It is found in most laurel forest plots, but 
also grows well at degraded sites and outside the range of 
the laurel forest ecosystems. M. faya is quite abundant due 
to its high colonizing capacity (Banares et al. , 2004). In ad- 
dition, Myrica faya has been recognized as one of the 
twelve most noxious plants alien to Hawaii due to its ability 
to rapidly and aggressively invade and colonize the Hawai- 
ian environment (Smith, 1985). 

M. faya co-occurs very close (separated by few me- 
ters) to Myrica rivas-martinezii A. Santos. The latter is an 
endangered endemic species of the laurel forests on the Ca- 
nary Islands. M. rivas-martinezii is found on rather poor 
soils along the inferior margins of the laurel forest domain, 
especially in isolated locations that have been the object of 
multiple exploitations. It is a perennial, woody species that 
was first time described in El Hierro in 1980 (Santos- 
Guerra, 1980). This dioecious tree species occurs on only 
three islands (El Hierro, La Gomera and La Palma) of the 
Canarian archipelago. Its most abundant population (with 
about 40 individuals) is located at El Fayal, in an area of ap- 
proximately 90 km 2 on El Hierro Island. Only ten isolated 
individuals are known in La Gomera Island, six of which 

Send correspondence to Miguel Angel Gonzalez-Perez. Departa- 
mento de Biologia, Universidad de Las Palmas de Gran Canaria, 
Campus Universitario de Tafira 35017, Las Palmas de Gran Cana- 
ria, Canary Islands, Spain. E-mail: mgonzalez@becarios.ulpgc.es. 



have male flowers, all in different locations and separated 
by linear distances ranging from 1 to 7 km, but isolated by 
deep ravines and cliffs. In La Palma, the only two known 
individuals (one male and one female) are separated by 
more than 20 km (Banares et al., 2004). Myrica rivas- 
martinezii has been classified as Critically Endangered ac- 
cording to IUCN categories (VVAA, 2000). It is also pro- 
tected by the Canarian Government (Boletin Oficial de 
Canarias - BOC, 2001) and the European Habitat Directive 
(Beltrane/a/., 1999). 

Previous studies showed no evidence for genetic dif- 
ferentiation on the basis of the sequences of either the 
pastid rbcl gene, which encodes the large subunit of rubis- 
co, or the 1 8S-26S nuclear internal transcribed spacer (ITS) 
region (Huguet et al, 2005). Yet, Batista et al. (2004) 
showed higher genetic diversity levels in M. rivas- 
martinezii using RAPDs. 

We sampled 41 specimens of Myrica rivas-martinezii 
and 183 of M. faya in El Hierro, La Gomera, and La Palma 
Island (Table 1). Microsatellite loci were isolated from M. 
faya genomic DNA using biotin/streptavidin enrichment 
(Khasa et al. , 2000). Briefly, genomic DNA was digested 
with restriction endonucleases (Haelll, Alul and Rsal) and 
linkers were added using T4 DNA ligase. Linker adapted 
total genomic fragments were then enriched by hybridiza- 
tion with 5' biotin AC 12 and TC12 followed by separation 
with streptavidin paramagnetic beads (M270S, Invitrogen, 
CA). After amplification with linker-specific primers, en- 
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Table 1 - Myrica rivas-martinezii and M. fay a populations analysed in the 
Canary Islands. 



Island 


Species 


N 


La Gomera 


M. rivas-martinezii 


10 




M. fay a 


132 


El Hierro 


M. rivas-martinezii 


30 




M. fay a 


38 


La Palma 


M. rivas-martinezii 


2 




M. fay a 


13 


Total 


M. rivas-martinezii 


42 




M. fay a 


183 



N = sample size. 

riched fragments were cloned into plasmid vectors 
(pGEM3Z+, Promega) and single colonies containing ACn 
or TCn microsatellites were identified by colony hybridisa- 
tion. After one round of biotin enrichment, 75/1000 ampi- 
cillin resistant colonies contained microsatellite sequences. 
Inserts from 23 positive colonies were amplified with Ml 3 
universal forward and reverse primers, treated with exo- 
nuclease I and Antarctic alkaline phosphatase (New Eng- 
land Biolabs) , and then sequenced from both orientations 
using a ABI3730 capillary electrophoresis system (NAPS 
Service, University of British Columbia). 

PCR primers complementary to the flanking regions 
of 13 loci were designed with approximately 40% GC and 
avoiding palindromic sequence motifs. Out of 13 micro- 
satellite loci analyzed, six yielded PCR products of ex- 
pected size and were considered robust and predictable 
enough for further analyses (Table 2). 

Total DNA was extracted following the Dellaporta et 
al. (1983) method modified by Corniquel and Mercier 
(1994). Each 25 uL PCR reaction contained approximately 
20 ng DNA, 10 pmol of each primer, 0.25 uL BSA (0.4%), 
as well as PCR Master Mix (Reddy-Mix, ABgene, Surrey, 
UK). Amplifications were carried out using the following 
conditions: 3 min denaturation at 95 °C, 35 cycles of 30 s 
denaturation at 95 °C, 30 s annealing at 55 °C, and 1.5 min 
elongation at 72 °C, followed by 5 min final elongation at 
72 °C. The products were detected using an ABI 3 100 Ge- 
netic Analyzer and fragment sizes were determined using 
GENESCAN v. 2.02 and GENOTYPER v. 1.1 (Applied 
Biosy stems, Inc.). We identified allele peak profiles at each 
locus and assigned a genotype to each individual. 

Exact Hardy-Weinberg tests to measure the signifi- 
cance of deviations from the null hypothesis of random un- 
ion of gametes (Guo and Thompson, 1 992) were carried out 
on natural populations of M. rivas-martinezii and M.faya 
using Fisher exact test, both for each pair of loci and within 
each population using GENEPOP 1.2. (Raymond and 
Rousset, 1995). Basic genetic diversity indices, mean num- 
ber of alleles (^4), the observed (Ho), and unbiased expected 
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(He) heterozygosities (Nei, 1978), paternity exclusion pro- 
bability (Q) (Weir, 1996), and probability of genetic iden- 
tity (T) (Paetkau et al, 1995), were estimated using 
IDENTITY 1.0. The combined probability of paternity ex- 
clusion, QC= 1 - [11(1 - Qi)] and the combined probability 
of genetic identity IC = ITY; were also estimated for overall 
loci. The proportion of null alleles for each locus was calcu- 
lated as (He-Ho)l\ +He, following Brookfield (1996). 

All six loci were polymorphic, the number of alleles 
ranging from 3 (M24) to 16 (M10), while Ho ranged from 
0.30 (M5) to 0.88 (M10), and gene diversity (He) from 0.32 
(M5) to 0.89 (M10). Only one locus (M24) showed values 
of proportion of null alleles over 10% in both species, sug- 
gesting the possibility of null alleles (Table 2). Except for 
locus M24, none of the other loci showed significant devia- 
tions from Hardy-Weinberg equilibrium in M. rivas- 
martinezii, (data not shown). The M24 locus was estimated 
to have a null allele frequency of 0.16. 

While there was no evidence for linkage disequilib- 
rium between loci in the data set of 225 individuals of both 
Myrica species, three (Mil, M20 and M24) of the six ana- 
lyzed loci showed significant defect of heterozygote in M. 
faya. These three loci were estimated to have null allele fre- 
quencies of 0.10, 0.04, and 0.17, respectively. Many of the 
loci showed high probability of paternity exclusion, which 
should provide powerful markers for paternity exclusion 
when genotype data from several loci are combined. Total 
paternity exclusion probability (QC) for this set of six 
microsatellite loci was estimated at 0.981 for M. faya, and 
0.941 forM. rivas-martinezii. Total probabilities of identity 
(IC) were 5.71 x 10" 6 and 9.03 x 10 -5 for M. faya and M. 
rivas-martinezii, respectively. 

Genetic diversity recorded for the endemic species M. 
rivas-martinezii was similar to those described for other en- 
dangered endemic species of the Canary Islands, for exam- 
ple Bencomia exstipulata (A = 6.20, Ho = 0.43, He = 0.65) 
(Gonzalez-Perez et al., 2004). In contrast genetic variabil- 
ity in M. faya was higher than corresponding values de- 
tected in other, not endangered native species of the Canary 
Islands, e.g., B. caudata (A = 6.40, Ho = 0.47, He = 0.62) 
(Gonzalez-Perez et al, 2004). This may reflect the high 
colonizing capacity of M.faya (Banares et al, 2004). 

Altogether, the primers described in this paper can 
provide useful markers to investigate genetic relationship 
between Myrica rivas-martinezii and M. faya. 
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